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Abstract 



1 Introduction 



The recent literature offers several results on the modeling of systems gov- 
erned by conservation laws on networks. For instance, in [5l El \12\ [T3] the 
modeling of a network of gas pipelines is considered. The basic model is 
the p-system or, in |21| . the full set of Euler equations. The key problem in 
these papers is the description of the evolution of fluid at a junction between 
two or more pipes. A different physical problem, leading to a similar ana- 
lytical framework, is that of the flow of water in open channels, considered 
for example in |33j . 

Recent papers deal with the control of smooth solutions, see for in- 
stance P3H2Zll2Hll2Sll3iai33lE5ll35ll35- 0ther approaches are based on 
suitable discretizations, as in [25l [391 ED E] • The present work presents a 
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general framework comprising several models in the existing literature and 
providing a proof of the existence of optimal controls for physically rea- 
sonable cost functions. In particular, in the structure below, the solution 
considered may well be non smooth and optimality is achieved in the set of 
all L 1 controls with bounded variation. 

As samples of the applications of the present results, we extend results 
from the current literature. First, we consider the optimal management of 
a compressor in a gas network. This device is required to compress fluid 
guaranteeing a given pressure while consuming the minimal energy, see \39\ 
I44j . Then, the present framework is used to cover different optimization 
problems for the flow in open canals: the keeping of a constant water level 
through the optimal management of an underflow gate, see [23J, and the 
prevention of overflow in a multiple valves system, see [32], or a pumping 
station, see [2*7] . 

From the analytical point of view, the above models are described by a 
system of conservation laws of the form 

t G [0,+oo[ 

dtui + d x fi(u{) = gi(t,x,ui) with x £ [0,+oo[ (1.1) 

I = 1, . . . , n. 

Here, u\ is the vector of the conserved variables along the Z-th pipe, fi is 
a general nonlinear flux function and g\ is the source term related to the 
'-th tube. A time dependent interaction at the junction is described by time 
dependent conditions on the traces of the unknown variables at the junction, 
namely 

(m(t, 0+), u 2 (t, 0+), . . . , Un(t, 0+)) = n(t) (1.2) 

for a suitable smooth ty, see [5j EJ [12] [13] [21] [32]. In the examples below, 
n(t) is the control to be chosen in order to minimize a given cost functional 
J = J{u,U). 

To obtain the existence of a control minimizing J ', we need first to prove 
the well posedness of (jl.l|) - ([1.2|) . This is the content of our main analytical 
result, namely Theorem 12 .31 The existence of an optimal control then follows 
in Proposition 12.41 

The next section is devoted to the analytical results. Section [3] presents 
the applications while all the technical details are collected in Section "3J 

2 The Cauchy Problem at an Intersection 

Throughout, we refer to [7J for the general theory of hyperbolic systems of 
conservation laws. Let SI/ C R 2 be a non empty open set. Fix flows fi 
such that / = (/i, . . . , /„) satisfies the following assumption at an n-tuple 
of states u = (ui, ■ ■ ■ , u n ) € 0, where SI = Q\ x 0,2 x • • • x Sl n : 
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(F) For I = 1, ... ,n, the flow ft is in C 4 (f^;R 2 ), Dft{u{) admits a strictly 

negative eigenvector X\(ui), a strictly positive one \ l 2 (ui) and each 
characteristic field is either genuinely nonlinear or linearly degenerate. 

Under this condition, when gi = and x G R, (jl.ip generates a Standard 
Riemann Semigroup, see [7J Chapter 8]. Recall that a 2 x 2 system of 
conservation laws admits entropies, see |43| Paragraph 9.3]. 

Here and in what follows, R + = [0, +oo[. For later use, with a slight 
abuse of notation, we denote by 

IMI = Sr=i ll^ll for u G Q 

||ti|| L i = J R+ 1 1 ti(ac) 1 1 dx for u G L^R+jfi) , 
TV(n) = £r=l TV W foru G BV(l+;0) . 

Below, the constant state u £ f!, is fixed. Throughout, we also fix a time 
T G ]0, +oo] and a positive 5. For all 8 G ] 0, 5] , we denote 

U s = {u G u + L X (R + ; 17): TV (u) < <$} . 

On the source term g = (gi, ■ ■ ■ ,g n ) we require that if G is defined by 
(G(t, u)) (x) = ^gi (t, x, «i(x)) , . . . , g n (t,x, u n (x)yj , then G satisfies: 

(G) G: [0,f] X i-> L^R+jR 2 ™) is such that there exist positive L 1} L 2 
and for all t G [0, f] 

\/u,w£Ug \\G(t, u) — G(s, iw)|| L i < £i • ( || ix — u>|| L i + |i — s|) 
VueWj TV (<?(*,«)) < L 2 . 

Below, we require only (G), thus comprising also non-local terms, see |18| . 
Examples of sources g such that the corresponding G satisfies (G) are pro- 
vided by the next proposition, which comprehends the applications below. 

Proposition 2.1 Assume that the map g:R + x $7 i— > R 2n satisfies: 

1. there exists a state u and a compact subset K, of R + such that g(x, u) = 
for all x G R+ \ K.; 

2. there exists a finite positive measure fj, such that for all x%,x 2 G R + 
with x\ < X2, and all ue!l', 

\\gi(x 2 +,u) -gi(xi-,u)\\ <^([x 1 ,x 2 }) ; 

3. there exists a positive L such that for all u, w G ft, for all x G R + , 

\\g(x, u) — g(x, w)\\ < L ■ \\u — w\\ . 
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Then, condition (G) is satisfied. 



The proof is deferred to Section 01 

We consider the Cauchy problem at a junction, see [HJ [I3j EUJ [22j [31] . 
First, we extend [20|, Definition 3.1] to the present case of a Cauchy problem 
with sources. 



Definition 2.2 Fix the maps * G C 1 (ft; R n ) and II G BV 

lueafc solution on [0, T] to 



dm + d x fi(ui) =gi(t,x,ui) 

*(«(t,o)) =n(t) 

it(0, x) = n (:c) 



f G R+ / G {1, . . . , re} 

x G R+ u G n + L 1 (R+;ft) 



(2.1) 



is a map u G C° ([0, T]; n + L 1 (M + ; f2)) sitc/t that for all t G [0,T], G 
BV(R+;0) and 

(W) u(0) = u and for all (p G C£° (]0, T[ x ]0, +oo[ ; R) and for I = 
1, . . . ,n 



n(ui d t <p + fi(ui) d x (p) dx dt+ I 
Jo 



<p(t, x) gi(t, x, ui) dxdt = . 



(ty) The condition at the junction is met: for a.e. t G R + , \£ (u(t, 0+)) = 

n(t). 

T/ie w/eafc solution (p,q) is an entropy solution if for any entropy - entropy 
flux pair (n h q{), for all ip G C£° (]0, T[ x ]0, +oo[ ;R + ) and for I = 1, . . . , n 



T 



il+ 



{r)i(ut) d t ip + qi (ui) d x ip) dx dt + / Dr)i(u{) g(t,x,u)tpdxdt >0 



o ii+ 



We are now ready to state the main result of this paper, namely the well 
posedness of the Cauchy Problem for (|2.ip at the junction. 

Below, we denote by r l 2 (u) the right eigenvector of Dfi(u) corresponding 
to the second characteristic family. 

As is usual in the context of initial boundary value problems, [H El [3], 
we consider the metric space X = (n + L 1 (R+,0)) x (tl + L 1 (R+,R n )) 

equipped with the L 1 distance. Let the extended variable p = (u, II) with 
u = u(x), respectively II = TL(t), be defined for x > 0, respectively t > 0. 
Correspondingly, denote 



d x ({u, n), (u, n)J = (u, n) - (u, n) = ||u - n|| L1 + 

TV (p) = TV (n) + TV (II) + II* (u(0+)) - 11(0+) 
V s = {p G V:TV(p) < 5} 

Below, 7j is the right translation, i.e. (7^11) (s) = U(t + s). 



n-n 



(2.2) 
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Theorem 2.3 Let n G N, n > 2 and assume that f satisfies (F) at u and 
G satisfies (G). Fix a map ^ G C^f^R") that satisfies 



det 



Di*(«)r^(«i) £> 2 ¥(w)r!(« 2 ) ... £> n ¥(«K(«») ^0 (2.3) 



where Di^ = D Ul ^> , and let II = ^(u). Then, there exist positive 5, 5' , L, T, 
domains T>t, for t G [0, T], and a map 

£:{(T,t ,p):t G [0,T[, r G [0, T - t ], p G V to } ^ P 5 

1. V s ' CPjCP s /or aM t G [0, T]; 
5. for all t G [0, T] and pGD t „, £"(0, i D )p = p; 
3. for all t G [0,T] and r G [0, T — t a ], S(r,t )V to C P to+T ; 
^. /or a// i D G [0, T] ; n, t 2 > with t\ + t 2 G [0, T — t ], 
£(r 2 ,t + ri) o £"(n,t ) = £(r 2 + n,t ) ; 

5. for all (u ,H) G X>t oJ set £(t,t )(u ,H) = (u(i),%U) where t *—>■ u(t) is 
the entropy solution to the Cauchy Problem ()2.1|) according to Defini- 
tion \2.S\ while the second component t i— > 7^11 is the right translation; 

6. £ is tangent to Euler polygonal, in the sense that for all t Q G [0, T], 
for all (n ,II) G V to , setting £(t,t )(u ,U) = (u(t),T t U), 



lim — 

t^o t 







i(t) - (St(u ,U)+tG(t ,u )) 
where S is the semigroup generated by the convective part in (|2.1[) : 
7. for all t a G [0, T], r G [0, T — t a ] and for all p, p G T> to , 



\£(t, t )p ~£{t, t )p\ 



< L-\\u 



+L 



to + T 



to 



nW-n(t) 



dt. 



(2.4) 



The proof is deferred to Section 0J Note that in the case n = 2, fx = — / 2 , 
II = and \&(ui,u 2 ) = + / 2 (n 2 ) we (re)obtain the well posedness of 

a standard 2x2 balance law. 

From 7. in Theorem 12.31 we immediately obtain the following existence 
result for an optimal control function II to the nonlinear constrained opti- 
mization problem 

( dm + d x fi(ui) =gi(t,x,ui) 
* (u(t,0)) = II(t) on [t a ,T] . 

u(0, x) = u {x) 



minimize J~(n) subject to < 
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Proposition 2.4 Let n G N, n > 2. Assume that f satisfies (F) at u and G 
satisfies (G). Fix a map ^ G C 1 (r2;R n ) satisfying ([23]) and /ei fl = ^(u). 
W^ii/i i/ie notation in Theorem \2.3\ for a fixed u Q (zUs, assume that 

Jo : |n |[oT] :n G (fl + L 1 (]R + ;]R n )) and (u D ,n) G P 5 | 
Jl : D* 5 i-» M 

are non negative and lower semicontinuous with respect to the L 1 norm. 
Then, the cost functional 

J(U) = J (n) + [ T Ji (£{t, 0)(u o , n)) dr (2.5) 
Jo 

admits a minimum on jll G (fl + L 1 ([0, T]; R n )) :(u ,II) G P }- 

Proof. Due to Theorem I2.3( f (-,0)u o is Lipschitz continuous in II. An ap- 
plication of Fatou's Lemma shows that the second summand in (|2.5p is lower 
semicontinuous. Hence, also J is lower semicontinuous and the existence of 
a minimizer follows from Weierstrafi Theorem. □ 



3 Networks of Gas Pipelines and Open Canals 

3.1 Compressor Control for Gas Networks 

We describe a compressor acting between two pipes with the same cross 
section by the equations 

d t pi + d x qi = 0, 

d t qi + d x \— +p(pi)) = -X rnT1 0) v^^- - gpi sin ai(x) 
\Pi J [°- L ] Pi 

t G M + , XGK+, I = 1, 2 , (pi,qi) G M + xR. 

where p is the mass density of a given fluid, q its linear momentum density, 
v accounts for the friction against the pipe's walls, g is gravity and a(x) is 
the inclination of the pipe at x. It is reasonable to assume that a{x) = for 
x sufficiently large. Furthermore, we are interested only in the dynamics in 
the pipe up to the maximum length L. The pressure law p = p(p) satisfies 

(P) p G C 2 (K+;R+), p(0) = and for all p G R+ p'{p) > 0, p"{p) > 0. 
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flow direction 



Tube 2 




Tube 1 



Figure 1: Notation for the compressor model (|3.1jl . 



As usual in the engineering literature, we focus below on the 7 law 



p(p) = V* 



(3.2) 



for suitable positive constants p*, p*. Choosing an initial datum u\ = (pi, qi) 
in the subsonic region 

U l = {(p,q) 6K+X R:\i(p,q) < < A 2 (p, ?)} , 

ensures that (F) holds at u. Recall the standard relations 

AiO,g) 



(?/p) - vV(p) > 



1 



-Ai(p,g) 



r 2(p,q) 



(q/p) + v¥{pj- 



1 



A2(p, 9) 



(3.3) 



The coupling condition describe the effect of a compressor sited at the 
junction between the 2 pipes. A standard relation in the engineering litera- 
ture is the following, see also [40, Section 4.4, Formula (4.9)] or 



#(iti,it 2 ) 



92 



qi + <?2 

P(Pl) 



1 



and U(t) 





n 2 (t) 



(3.4) 



and II2 is proportional to the applied compressor power. 
Proposition 3.1 Let «i,«2 £ ^ satisfy 



for a positive 1I2 . Then, Theorem \2.3\ applies to (|3.ip — (|3.2[) -([3T 





n 2 



Proof. By (I3.3p . (F) is satisfied at it. Condition (G) is satisfied due to 
Proposition 12.11 Finally, condition (|2.3j) leads to the determinant 



det 





1 




1 






d pi ^2 









1 

9 P2 *2 9«j 2 *2 



1 

A2(« 2 ) 



A 2 («i) (A 2 (^ 2 )3 52 *2 + d P2 V 2 ) - X2{u 2 )d Pl ^ 2 > 



7 



due to the choice ui,U2 £ Hence, Theorem 12.31 applies. 



□ 



A typical optimization problem in gas networks [391 l4"4"] is the control 
of compressors stations such that a certain outlet pressure p for a customer 
located in the interval [x a , x&], x a > 0, is satisfied. In the optimization 
problem we penalize large energy consumption by the L°°-norm, frequent 
changes in the applied compressor energy by the TV —norm and deviations 
from the desired outlet pressure. We model this situation by considering 



j (ii) = TV(n) + ||n| 

Jt (£(r,0)(« o ,II)) = 



x b 



and 

\p(P2(t,x)) -p\dx, 



where p2(t,x) is given by the solution of (13. ip . The lower semicontinuity of 
J is obvious. J\ is L 1 -Lipschitz. Indeed, let p2 and p~2 denote the density 
distribution in the pipe I = 2 corresponding to the same initial datum and 
to the controls II and II. Then, 



Jt (£(r,0)K,II)) - Ji (S(t,0)(u o ,R) 

T f -x b 

|p(P2(*o + t, x)) -p\- \p(p~2(to + r, x)) - p\ 



< 



JXa 

rT rx b 

<C I I \p2(t, x) — p~2(r,x)\ dxdt 



dx dr 



< C 



J X a 

n 2 - n. 



Li([0,T]) 

for some constant C and due to (|2.4p . Hence, Proposition 12.41 applies. 



3.2 Control of Open Canals 

Similarly to the models in [23 |. 127 ^ 133]. we consider canals with fixed rectan- 
gular cross section having width b\ described by 



d t H t + d x Q t = 



n , , Qi\Qi\ ( 3 - 5 ) 

-gHismai(x) - X[0,h]i x ) v H 



where 



t 6 



x€R + , l = l,...,n, (H h Qi)€ 



and Hi(t,x) is the level of water at time t, point x in canal biQi is the 
total water flow; a\ is the inclination, g gravity and L is the length of the 
canal as in Section [3TT1 Note that (|3.ip reduces to (|3.5p in the case 7 = 2. 
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3.2.1 The Case of an Underflow Gate 

Following [23], consider (13, 5p for n = 2 with the coupling condition 



*(ni,n 2 ) 



£>iQi + hQ2 



and U(t) 





u{t) 



(3.6) 



Hi-if 2 

the control u being the opening of the underflow gate, see Figure EJ The 



Figure 2: Notation for (|3.5p with coupling conditions (|3.6|) . 
conditions (F) and (G) are proved as above. The determinant in (|2.3p gives: 



det 



6i 

2Qi 

{H1—H2) 2 H!-H 2 





Q1 2 



1 

A 2 (*ui) 



b 2 




1 

A 2 (w 2 ) 



(6iA 2 («i) + 6 2 A 2 (n 2 ))Qi 2 26 2 A 2 (ni)A 2 (u2)Qi 



(#1 - ff 2 ) 



> 0. 



The determinant is positive, since the underflow gates are only operating for 
Hi > H2 and (due to our parametrization of the pipe) for Qi < 0. Finally, 
Theorem 12.31 applies for sub-critical data u. 

A typical problem for horizontal pipes is to maintain a steady height in 
the downstream canal I = 2. Hence, we consider equation (13. with ct\ = 
and penalize large gradients in the water height. We introduce the cost 
functionals 

i-T p+00 

J = \u(t)\dt and J x = / <p(x) d\d x H 2 \ . (3.7) 

Here, ff 2 is the water height in canal I = 2 given by the solution to (|3.5p 
and (|3.6j) . The non negative and lower semicontinuous weight ip assigns 
different importance to oscillations in the water level at different locations. 
Under the assumptions of Theorem 12.31 the map x — > // 2 (t, x) is a function 
of bounded variation and hence, d x H2(t,x) is a Radon measure. Then, 
the measure \d x H2\ is the total variation of d x Hi- Due to (|2.3|) we ob- 
tain that n& — ► n* in L 1 implies that for any fixed t € [0, T] we have 
flf(V) H^(t,-) in L 1 . Therefore, the same arguments as in |14[ The- 
orem 2.2] and [144 Lemma 2.1] show that Proposition 12.41 can be applied 
to (|33 ]) -(|57f ]) . 
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3.2.2 The Case of Multiple Valves 

We consider conditions for valve control similar to those introduced in [9j 
Section 2.9] or [10] and discuss a situation with n connected pipes as in 
Figure We control the inflow at each connected pipe by the opening of a 



flow control valve 



This amounts to 



E?=i b iQi 

W(ui, . . . ,u, 

Qn-l 

The assertions of Theorem 



and n(i) 





ui{t) 



(3.8) 



u n -i(t) 

are satisfied, since the determinant in (|2.3p 




Figure 3: Illustration of a multiple valves junction. The outflow in the large 
canal is controlled through n — 1 valves at the incoming pipes. 

n 

evaluates to \\ bi\2(ui) ^ for any sub-critical state u. 
i=l 

We consider the problem to prevent overflow in the downstream canal n 
by valve control at the node. We assume costs associated with the operation 
of valves given by non-negative, bounded functions Cj(i) for i = 1, . . . , n — 1 
and a maximal height of the water h. We model this problem by minimizing 



n-1 .t r-l 

J = / Ci(t)ui(t)dt and J\ = j 
Jo Jo 



h) dx. 



(3.9) 



Herein, H n is the solution to (|3.5|) and (|3.8|) and x + = max{0,a;}. Since 
x — > x + is Lipschitz, we can apply the same arguments as in Section \3. II and 
obtain that Proposition 12.41 applies in this case. 
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3.2.3 The Case of a Pumping Station 



We consider the case of a simple pumping station, i.e. we supply (13. 5p with 
n = 2 and with the coupling conditions from 



*(«i,it 2 ) 



Hi — H2 



U(t) 





u(t) 



In the present case, conditions (F) and (G) are proved as in Proposi- 
tion [3TTJ A direct computation allows to verify that the determinant in (|2.3[) 
is non-zero and therefore Theorem 12.31 applies. A reasonable optimization 
problem consists in minimizing J for (j3.9f) . 



4 Technical Details 



As a general reference on the theory of hyperbolic systems of conservation 
laws, we refer to [7]. As usual, C denotes a sufficiently large constant de- 
pendent only on / restricted to a neighborhood of the initial states. 

Proof of Proposition \2.H Note first that 1. and 3. imply that G attains 
values in L 1 . Concerning the bound on the total variation, by 2., 

TV(Cr(u)) = sup ^ ff (xj, u(xj)) - g (a?j-i, u(xj-x)) 

i 

< sup ^2 5 ( X *> U ( X i)) ~ 9 (xi-l,u(Xi)) 

i 

+ sup^ L (xi~i,u(xi)) - g (xi-\,u(xi)) 

i 

< sup'Y]iJ,([xi-i,Xi]) + L sup^ \\u(xi) - u(Xi-l)\\ 

i i 

< 2/i(R + ) + LTV(ii). 

Condition 4. directly implies the L^'-Lipschitz condition on G. □ 



4.1 The Convective Part 

This section is devoted to the Cauchy problem for fll.ip — (|1.2|) in the case 
gi = for I = 1, ... ,n. First, we rewrite it as an initial - boundary value 
problem for a (2n) x (2n) system of hyperbolic conservation laws. To this 
aim, introduce positive 



A mltl < min inf 



and A' 



> max sup 

i=1 > 2 «efii 



Al(«) 
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for I = 1, . . . , n. Introduce the flow 

F 2l _ l {U) = A l -{f l ) l {U 2l ^U 2l ), F 2l {U) = A r {fi) 2 (U 2 i-i, U 2l ) (4.1) 
where the dilatation factor A/ are recursively defined by 
1 



Ai 



min 



and A 



A/_i for I = 2, . . . ,n 



Lemma 4.1 The flow (|4.ip defines a hyperbolic (2n) x (2n) system of con- 
servation laws, the eigenvalues Ai, . . . , A 2n of DF satisfy, for I = 1, . . . , n, 



A, 1 



1 \ k 

n ^max 

k=l mm 

n A max 
A fc . 

fc=l mm 



" < Ai < 
i 

1 < A 2 < 

< A 2i _i < 

A 2 z 



-1 



min 



< 



< 



n ^max 
A fc . 

fe=l mm 

' \ k 

n A max 
\k 

k=l mm 



Proposition 4.2 T/ie Cauchy problem at the junction and the IBVP 



dm + d x fi(ui) = 

tt(u(o+,t)) = n(t) 

w(0, cc) = u (x) 



6oi/i defined for t > and x > 0, mt/i 
(^ok-iM = (lt )2Z-l(x/Az) 



a t c/ + a^c/) = o 
*(tf(o+,t)) =n(t) 

U(0,x) = U o {x) 



(U ) 2 i(x) = (u ) 2l (x/Ai) 



are equivalent, in the sense that u = u(t, x) solves the former problem in the 
sense of Definition \2.S\ if and only if the map U = U(t, x) defined by 

U 2 i-i(x) = u 2 i-i(x/Ai) U 2 i{x) = u 2i (x/Ai) 

is a weak entropy solution to the latter problem. 

For the definition of weak entropy solution to the IBVP above, see }26] 
or [TJ El [2l]. The proof of Proposition 14.21 is immediate. Note that the 
two problems differ by a linear change of coordinates in the space variables, 
hence the entropicity of solutions is maintained. 
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4.1.1 The Riemann Problem at a Junction 

Let / satisfy (F) at u and let G C 1 (r2;IR n ). By Riemann Problem at the 
Junction we mean the problem 

dm + d x fi(ui) = , , 

■uz(0,x) = u ,z , 

where, for I = 1, . . . , n, u ,i are constant in Q l and II G M. n is also a constant. 

Definition 4.3 Fix the map ^ G C 1 (Q;IR Tl ). ^4 solution to the Riemann 
Problem ()4,2p is a function u: R + x IR + i— > SI sitc/i i/iai 

(L) For I = 1, ... , n, £/ie function (t,x) *— > ui(t,x) is self-similar and coin- 
cides with the restriction to x > o/ i/ie Lax solution to the standard 
Riemann Problem 



dm + d x fi{ui) = 

u 0j z ifx>0 
«i(l,0+) ifx<0. 



ui(0,x) 



(ty) The trace u(t,0+) of u at the junction satisfies (u(t, 0+)) = II /or 
aZZ f > 0. 

The following proposition yields the continuous dependence of the solution 
to the Riemann problem from the initial state, from the coupling condition 
and from the control term IT. 

Proposition 4.4 Let n G N with n > 2, and f satisfy (F) at u. Fix 
* G C 1 (r2 n ;M n ) satisfying (J2~3l) and a constant p = (u,TV) with U = ^f(u). 
Then, there exist positive 5, K such that 



1. for all p = (u ,H) with ||p — p|| < 5, the Riemann Problem (|4.2p 
admits a unique self-similar solution (t, x) i— » (TZ(p)) (i, x) in the sense 
of Definition ^. 2\ 

2. let p, p satisfy ||p — p|| < 5 and ||p — p|| < 5. Then, the traces at the 
junction of the corresponding solutions to (|4,2p satisfy 

(K(p)) (t, 0+) - (K(p)) (t, 0+)|| <K-\\p- p|| ; (4.3) 

3. call S(p) the n-vector of the total sizes of the 2-waves in the solution 
to Then, 

||E(p)-S(p)|| < if • ||p - p|| . 

The proof is omitted, since it follows from Lemma 2.2] through Proposi- 
tion 22] or from simple modifications of [20, Proposition 2.2]. 
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4.1.2 The Cauchy Problem at a Junction 

For a piecewise constant function u = u a v a _ 1 the usual Glimm 

functionals in the case of a non characteristic boundary, see [21 Lemma 4] 
and [U [M] , take the form 



V(u) 

Q{u) 
T(p) 



E 
E 



2Kj- 



a 



l.tt 



+ 



»,£» j,/3 



V(u) + i^T • TV (II) + K ■ Q(u) 



(4.4) 



where A 1 denotes the set of approaching waves in the Z-th pipe, see [7J 
Paragraph 7.3], while a\ a is the (total) size of the i-wave in the solution of 
the Riemann problem at x a in the Z-th pipe. Note that at x a = 0, we consider 
the Riemann problem at the boundary, according to Section 14.1.11 The 
constant Kj is defined as in |20| formula (6.2)], K is as in [20, Paragraph 6] 
and K is as in [21 Section 6]. 

The lower semicontinuous extension of T to all functions with small total 
variation is achieved in [19] in the case of the Cauchy problem on the whole 
real line. Here, we use the analogous result on the half line x > and the 
lower semicontinuity of the total variation with respect to the L 1 norm, see 
[15] for details. 

Moreover, in the proof of the Lipschitz continuous dependence of (|4.7|) 
with respect to the initial datum and the condition at the junction, an 
excellent tool is the stability functional introduced in [HI EH [3H] , see also [71 
ES]: 



*(p>p) = £E 



i=l 1=1 



r+oo 








Jo 





n-n 



(4.5) 



where u, u are piecewise constant functions in U$ and we let (q[(x), q l 2 (x) 
q (u l (x),u l (x)), see Chapter 8]). The weights W\ are defined by 



W\(x) 
W\{x) 



K-il + m A[(x) + K\ k 2 (T(p) + t( p ; 

1 + Kl A l 2 ( X ) + Kt K 2 (T(p) + T(p)) 



(4.6) 



for suitable positive constants k\ , k 2 defined similarly to [71 Chapter 8] , see 
also [19] . and K as in [201 Section 6]. Here, T is the functional defined 
in (|4.4p . while the A\ are defined by 



A\(x) 
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x a < x, a G i7(w) 
x Q > x, a G i7(w) 

x a < x, a G 
x a > x, a G 



if q\{x) < 0, 
if ql(x) > 0, 



see [71 Chapter 8] . Note that the lower semicontinuous extension of 3> to all 
functions with small total variation, defined in [TU], keeps all the properties 
of the original functional In particular, there exists a constant C such 
that for all p, p G T>$, 



1 

C 



\U - U T i + 



n-n 



J <*(p,p)<C(||«-fi|| L i 



n-n 



Proposition 4.5 Let n G N, n > 2 and / satisfy (F) at u. Let YL = ty(u). 
Then, there exist positive 5,L and a semigroup P: [0, +oo[ x V i— > D snc/i 



5. 



1. T> ~2 cl L iP d ; 

& /or a// (u, n) G V, P t (u,U) = (S t {u,U),T t U) , with P p = p and for 
8,t>0, P s P t p = P s+ tP; 



3. for all (u ,n) G T>, the map t \— > St(tt ,n) solves 

< 



<?t« + d x f(u) = 

*(«(o+,t)) =n(t) 

u(0, x) = "u (x) 



(4.7) 



according to Definition s. 1 ^ 
4- for p,p£P and t, i > 



l^tp — <S<p|| L i( 
|5ip-5 t -p|L lf 



-) 



-) 



< L 



< L 



\ u ~ ^IIl 1 (r+) + 



n-n 



Li([0,t]) 



i - 1 



5. if p G D is piecewise constant then, for t > sufficiently small, StP 
coincides with the juxtaposition of the solutions to Riemann Problems 
centered at the points of jumps or at the junction; 

6. for all p G T>, the map t \— > T(Pjp) is non increasing; 

7. for all p,p£D, the map t i— > 3>(-PtPi Pp) * s n on increasing; 
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there exist constants C, rj > such that for all t > 0, for all p, p G T> 
and v G L 1 (M+; ft) wtfi TV (v) < f] 



\S t p - S t p - v\\ L i 



< L I \\u — u — v\\t 1/ 



+ 



n-n 



Li([0,i]) 



+ C ■ t ■ TV (v) . 



Thanks to Proposition 14,21 the above result falls within the scope of the 
theory of initial - boundary value problem for hyperbolic conservation laws, 
see [H El [3] and p2] for the details. 

We now consider the source term. 

Proposition 4.6 Let g satisfy (G). For all (n, II) G V, for all t Q G [0, T] 

and all r > sufficiently small, the following relation holds: 

T(u + rG{t , «),n) < T(n, II) + C ■ r 
* + rG(t , u),IL),(u + rG(t , u), Ujj < e Cr & ((«, n), (fi, fi)) 

Proof. We consider piecewise constant functions, leaving the lower semi- 
continuous extension to general functions to the techniques used in the proof 
of pH Lemma 2.3]. 

For piecewise constant functions, the proof is obtained through a careful 
control of interactions. Off from the boundary, the same computations of |18} 
formula (3.6), Lemma 3.6] or [H Lemma 2.2] hold. 

Due to the junction, we have one more term: 

T (u + rG(t ,u),U) -T(«,n) 
= V (u + rG{t , u),U) - V(u, n) + K (Q (u + rG(t , u),U) - Q{u, n)) 



< Ct + Ct 

< Ct 



{G(t ,u)) (0+) 



The first estimate is thus proved. Concerning the second one, note that the 
q\(x) are not affected by the presence of the junction, so that exactly the 
same computations in the case of [TTJ, [18] hold. □ 

With the constant C defined above, let 

P* = {(u 5 II) eP:T((«,n)) <8-C(T-t)} 
where T <T and T < 6/C. 
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We now aim at showing that the present situation falls within the scope 
of [16]. In the metric space X equipped with the L 1 distance ()2,2p introduce 
the local flow 



F Mo (u,II) = (S t (u,Tl)+tG(t ,St(u,Tl)) 
In the next proposition, we refer to Definition 2.1 and Condition (D)]. 
Proposition 4.7 F is a local flow. 

Proof. For t Q G [0, T], r G [0, T - t Q ] sufficiently small and u G V to , due to 
the first inequality in Proposition 14.61 we have 

T(F Tjio (n,n)) = T^S T (u,U)+rG(t ,S T (u,U)),T T U(to)) 

< T(P T (n,n)) +Ct 

< T(u,U) + Ct 



and therefore 



F T , to {V to ) QVta+T- 



We now prove the Lipschitz dependence of F from r and (it, II). Note 
that (G) implies the boundedness of G for t Q G [0, T] and u G T>t a . Let 
t~i,T2 G [0, T — t Q ], then, by 4. in Proposition 14.51 and (G) 

d x (F Tuto ( u ,n),F T2 , to (u,n)) < 
< || s n (u, n) - 5^(14,11)11^ 



+ ||T ri n-T T2 n|| L1 



+ t\G (t , s T1 (u, n)) - t 2 g (t , s T2 (u, n)) 

< C-|r 2 -ri| + CTV(n)|T 2 -ri| 

< C-\t 2 -tx\. 



Similarly, we get for t Q G [0, T], r G [0, T — t Q ] sufficiently small and 

(u,n),(£t,n) g v 



dx [F T>to (u,U),F Tito (u,U] 



+ 



< S T (u,U) -S T (u,U 
T T U - T T fl 

< C ■ \\u — tt|| L i + 
completing the proof. 



+ T 



G(t ,S T (u,U)) -G(to,S T (u,U) 



L 1 



n-n 



□ 
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Lemma 4.8 There exists a constant C such that for all t Q G [0, T], for all 
(u, II), (u,II) G V and e > sufficiently small 

*(F £)to (n,n),F Mo («,fi)) <(i + ce)*((«,n),(w,n)) . 

This proof follows directly from Proposition 14.61 and 7. in Proposition 14.51 

For e > 0, recall the definition of the Euler e-polygonal F £ generated by 
F, see [HH Definition 2.2]. Let = [r/e], [•] denoting the integer part. 



k— 1 
h=0 



(4.8) 



The conditions that allow to construct a process generated by F that yields 
solutions to are proved to hold in the following proposition. 

Proposition 4.9 The local flow F satisfies the conditions: 

1. there exists a positive C such that for all t Q G [0, T], r € [0,T — t ], 
(u, IT) G V to and all k G N w#i (A; + l)r G [0, T - Q, 

dx (F h)io+T oF Ti ,>,n),F (H1)Tii J«,n)) <CkTT- (4.9) 

2. i/iere exists a positive constant L such that for all e G [0,6], for all 
t Q G [0, T\, t > sufficiently small and for all (it, II); (u, II) G Pt D 



d x ^ o (n,n),F^ o (n,n)j <L-d x ^,II),(fi,n)J . 
Proof. Consider the two conditions separately. 

1. (We follow here the same steps in the proof of |17} Theorem 1.1]). We 
refer to the map G introduced in (G). G is L^bounded, L 1 -Lipschitz and 
TV (G(u)) is uniformly bounded for u G V. Moreover, using 4. and 8. in 
Proposition 14.51 



x 



X 



F kT j 0+T F Tito (u,lT) - F( fc+1)T>io («,II) 
Fkr,t +r {S T P + rG(t Q , S T p),T T U) 

- (%+i)rP + (k + l)rG{t , 5( fc+ i) r p), T (fc+1)r n 
Skr {S T P + rG(t a , S T p),T T U) 
+krG (t + t, S kT (S T p + rG(t , S T p),T T U 
%+i)rP + (k + l)rG(t , S {k+1)T p 



+ 



T kT T T H — 7( fc+1 ) T II 



L 1 
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< 



< 



Skr (SVp + rG(t , S T p),T T U) - S kT (S T p, T T U) - rG(t , S( k+1]r p) 



+kr 



G[t + T, S kT (S T p + rG(t , S T p),T T U) - G(t a , S {k+1)rP 



L 1 



L 1 



rG(t , S T p) - rG(t , S {k+1)T p) L1 +Ckr TV (rG(t , S T p)) 



+krL 1 t + 



S kr (S T p + rG(t , S T p),T T U) - S kT (S T p, T T U) ^ 



x 



< Ct\\S t p — S kT (S T p, T T H)\\ L1 + C L2 kr r 

+kr L\ ■ (T + C\\rG{t ,S T p)\\^ 

< Ckrr. 

2. Let k = [r/e] with k € N. By Lemma 14.81 and 7. in Proposition 14.5 
Fl to (u,U) - Fl to (u,ti) 

< C*(F^ o (n,n),F^ o (fi,n) 

< c (l + C(r-ke)) *(F fc £ £)to («,n),F^ o (u,n 

< C (i + Ce) [i + c(T-ke)) *(i^_ 1Kto («,n),j^_ 1)eito (fi,n) 

< ... 

< Ce Cr * f(u,n),(w,n)' 



< Ce 



CT 



(n,n) - («,n) 



completing the proof. 



□ 



Proof of Theorem \2.tA The application of [TBI Theorem 2.5] yields a 
process £, Lipschitz with respect to the L 1 distance on X, proving 1. to 4. 
in Theorem 12 .31 

The tangency condition 6. follows from [161 (2.9) of Theorem 2.5], indeed 



L 1 



u(t) - (S t (u ,U) + tG(t ,u )) 
< \\St,t {.u ,U) - F tt t (u ,Il)\\ x 

F tito {u ,U) - (S t {u ,U)+tG{t ,u ),T t IL) 



x 



< Ct 2 + t 

< Ct 2 . 



G(t ,St(u ,IL))-G(t ,u ) 



L 1 



Conditions (W) and are an easy consequence of the tangency condition, 
see pS Theorem 1.2]. 
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Finally, by [161 b) in Theorem 2.5], the process £ is Lipschitz, i.e. 

|f(r,t )(«,n) - £(T,t )(u,U) 

The better bound (12. 4h follows from the above construction. Indeed, the 
approximate solution F^ t (u, II) depends only on the restriction of II to 

[t ,t + T]. □ 
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